hcp-ph/9709281 



PITHA 97/35 
1 



Next-to-Leading Order QCD corrections 
to Three- Jet Rates with Massive Quarks* 

W. Bernreuther^ , A. Brandenburg^, and P. Uwer 



Oh: 

00 



00 

o 

Oh: 



X 



Institut f. Theoretische Physik der RWTH Aachen, D-52056 Aachen, Germany 

The reaction e'^e~ annihilation into three jets was recently computed for massive quarks at next-to-leading 
order in perturbative QCD. We discuss some results of our calculation for b jets produced at the Z resonance. 



1. INTROCUCTION 

Huge samples of Z decays into jets of hadrons 
have been collected both at LEP and at SLC. 
These data allow for numerous precision tests of 
the theories of electroweak and strong interac- 
tions. In particular, large numbers of jet events 
involving b quarks can be isolated using vertex 
detectors. For detailed investigations of b jets, 
aiming at a precision at the few percent level, one 
must take into account in theoretical predictions 
that the b quark mass is non-zero. In fact b quark 
mass effects in jets produced at the Z resonance 
have already been seen indirectly in the so-called 
"flavour-independence tests" of the strong cou- 
pling. Such tests have been performed both at 
LEP § and at SLC §. 

Three groups have recently reported [^-^ the 
calculation of the next-to-leading order (NLO) 
QCD corrections for e^e~ —> 3 jets for massive 
quarks. This extends the well-known NLO 3-jet 
results for massless quarks obtained in the early 
eighties |^^. (For the NLO 2-jet cross section 
and leading order (LO) results for three, four, and 
five jet rates involving massive quarks see and 
10, respectively.) In this talk we discuss some 
results of our NLO computation of (Tsjot and ap- 
ply it to b jets produced at the Z resonance. 
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2. THE NLO CROSS SECTION 

The calculation of the e^e^ annihilation cross 
section, <j%i^qj into three jets involving a massive 
quark antiquark pair to order consists of two 
parts: First, the computation of the amplitude 



of the partonic reaction e'^e 
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QQg 



at leading and next-to-leading order in the QCD 
coupling. Here Q denotes a massive quark and 
g a gluon. We have calculated the complete de- 
cay amplitude and decay distribution structure 
for this reaction. This allows for predictions in- 
cluding oriented three-jet events. The differen- 
tial cross section involves the so-called hadronic 
tensor which contains five parity-conserving and 
four parity-violating Lorentz structures. Second, 
the leading order matrix elements of the four- 
parton production processes e'^e^ Z*,^* 
ggQQ,QQqq,QQQQ are needed. Here q denote 
light quarks which are taken to be massless. 

The infrared (IR) and ultraviolet (UV) singu- 
larities, which are encountered in the computa- 
tion of the one- loop integrals, are treated within 
the framework of dimensional regularization in 
Z? = 4 — 2e space-time dimensions. We remove 
the UV singularities by the standard MS renor- 
malization. We have converted from the outset 
the on-shell mass of the heavy quark Q into the 
corresponding running MS mass. It is known that 
far from threshold one thereby absorbs some large 
logarithms into the running mass. 

An essential aspect of any NLO computation 
of jet rates is to show that the IR singularities of 
the virtual corrections are cancelled by the singu- 
larities resulting from phase space integration of 
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the squared tree amplitudes for the production of 
four partons. Different methods to perform this 
cancellation have been developed (see [p"2|-|l^ and 
references therein). We use the so-called phase 
space slicing method elaborated in The ba- 
sic idea is to "slice" the phase space of the four 
parton final state by introducing an unphysical 
parton resolution parameter Smin <C sycut, where 
Ucut is the jet resolution parameter. The parame- 
ter Smin splits the phase space into a region where 
all four partons are "resolved" and a region where 
at least one parton remains unresolved. For mass- 
less partons, the resolved region may be conve- 
niently defined by the requirement that all invari- 
ants Sij = {ki + kjY constructed from the parton 
momenta ki are larger than the parameter Smin- 
We have modified this definition slightly to ac- 
count for masses. 

In the unresolved region soft and coUinear di- 
vergences reside, which have to be isolated ex- 
plicitly to cancel the singularities of the virtual 
corrections. This is considerably simplified due 
to coUinear and soft factorizations of the matrix 
elements which hold in the limit Smm — > 0. (In 
the presence of massive quarks, the structure of 
coUinear and soft poles is completely different as 
compared to the massless case.) After having 
cancelled these IR poles against the IR poles of 
the one-loop integrals entering the virtual correc- 
tions, one is left with a completely regular differ- 
ential three-parton cross section which depends 

The contribution to cr^n^o "resolved" 
part of the four-parton cross section is finite and 
may be evaluated in Z) = 4 dimensions, which 
is of great practical importance. It also depends 
on Smin and is most conveniently obtained by a 
numerical integration. Since the parameter Smin 
is completely arbitrary, the sum of all contribu- 
tions to CTpfj^Q^ycut) must not depend on Smin- In 
the soft and coUinear approximations one neglects 
terms which vanish as Smin — *■ 0. This limit can 
be carried out numerically. Since the individual 
contributions depend logarithmically on Smim it 
is a nontrivial test of the calculation to demon- 
strate that (Jpfj^Q becomes independent of Smm 
for small values of this parameter . 

The three jet cross section depends on the 



experimental jet definition. We consider here 
the JADE and Durham clustering algo- 
rithms, although other schemes can also be 
easily implemented. We have checked that we re- 
cover the result of in the massless limit. 

Let us now discuss some results for the cross 
section ct^^q for b quarks. An important point to 
notice is that we are concerned here with the com- 
putation of tagged cross sections - whereas the 
massless NLO results [^j-^ apply to summing over 
all quark flavours in the final state. In the com- 
putation of tagged cross sections one encounters 
mass singularities in the real corrections - that 
are regulated by keeping the quark mass non-zero 
- which find no counterpart in the virtual correc- 
tions against which they can cancel. In the case at 
hand, to order a^, there is a contribution to ct^^q 
from the diagrams e^e~ —^ qqg* —^ qqbb. This 
contribution contains large logarithms involving 
irib- (They result from the region where the in- 
variant mass of the virtual gluon becomes small.) 
One may simply keep these logarithms, but that 
is not satisfactory, in particular if ^/s 3> mi,. 
Eventually, one has to factorize these logarithms 
into a fragmentation function for a gluon into a B 
hadron. A detailed discussion of the single 6-tag 
cross section ct^^q will be given elsewhere |l9| . 

Here we shall consider instead the following 
"double &-tag" three-jet cross section: we require 
that at least two of the jets that remain after 
the clustering procedure contain a 5 or 6 quark. 
The cross section ctji^^q, where the contributions 
from the above diagrams are included, remains 
infrared-safe in the limit mf, — > 0. 

As mentioned above we have expressed ct^^q 
in terms of the b quark mass parameter mh{ii) 
defined in the MS scheme at a scale ^. The 
asymptotic freedom property of QCD predicts 
that this mass parameter decreases when be- 
ing evaluated at a higher scale. (A number 
of low energy determinations of the b quark 
mass have been made; see for instance pO|-|2^ 
and references therein.) With mb[mf,) = 4.36 
GeV ll) and a,{mz) = 0.118 ||2|l as an in- 
put and employing the standard rcnormalization 
group evolution of the coupling and the quark 
masses, we use the value mf,(/i = mz) = 3 
GeV. These values for as and mb are used in 
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Figs. la,b where we plot ct^^q as a function of 
Vcut together with the LO result at the Z peak, 
both for the JADE and the Durham algorithm. 



by means of differential two-jet distributions, for 
instance in the Durham scheme. 





Figure 1: The cross section a^''' as a function of ycut for the JADE and Durham algorithms, respectively. 
The dashed line is the LO result. The NLO results are for renormalization scale ji = mz (solid line), 
fj, = mz /2 (dotted line), and /i = 2m z (dash-dotted line). 



(Initial state photon radiation and electroweak 
corrections are not included.) The QCD correc- 
tions to the LO result are quite sizable, as is 
known also in the massless case. At small val- 
ues of 2/cut) for instance below ycut ~ 0.01 for 
the JADE algorithm, perturbation theory is not 
applicable. The NLO corrections reduce the de- 
pendence on the renormalization scale significi- 
natly, as compared to the leading oder result for 
the three-jet cross section. This renormalization 
scale dependence, which is also shown in Figs. 
la,b, is modest in the whole ycut range exhibited 
for the Durham and above ycut ^ 0.01 for the 
JADE algorithm. 

The effects of the non-zero b quark mass on 
the three-jet rate at the Z peak is of the order 
of a few percent, depending on the jet algorithm 
and on ycut- An interesting application is the de- 
termination of the MS mass parameter at a 
high scale by measuring a suitable double ratio of 
three-jet fractions, as proposed in ||2^, computed 
to NLO in 1^ , and experimentally pursued by the 
DELPHI collaboration In we have com- 
puted a double ratio that is somewhat differently 
defined than the one used in Effects due 

to the non-zero b quark mass may also be studied 



3. CONCLUSIONS 

We have computed the complete differential 
distributions for e^e~ annihilation into three and 
four partons, including the full quark mass depen- 
dence, both on and off the Z resonance. Apart 
from cr^5.o tti^se results allow the computation of 
a number of other observables at order a^. Our 
results may be applied to b and c quark jet pro- 
duction at various cm. energies, and to theoreti- 
cal investigations of top quark production at very 
high-energetic electron positron collisions. 
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